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Using the example of graphene, we have extended the classic ripple tank experiment to illustrate the behavior 
of waves in periodic lattices. A loudspeaker driving air through a periodically perforated plexiglass plate 
onto the tank’s water surface creates wave patterns that are in agreement with numerical simulations and 
are explained in terms of solid state theory. From an educational point of view, the experiment provides an 
illustrative example of the concepts of reciprocal space and symmetry. 


At public outreach events as well as in schools or un¬ 
dergraduate lectures, the ripple tank is widely used to 
illustrate wave phenomena.I^^ However, a standard rip¬ 
ple tank’s one or two dappers cannot create patterns that 
illustrate the behavior of waves in periodic lattices. This 
challenged us to modify the setup in a way that provides 
multiple periodic sources. Employing a loudspeaker and 
a plexiglass plate with hexagonally arranged holes we 
were able to realize wave patterns in the ripple tank, 
that are reminiscent of the microscopic images of artif- 
ical graphene.l^ Two characteristic interference pattens 
are depicted in Figs, [^a) and (b). Readers specialized 
in solid state physics will recognize the same symmetry 
as graphene electronic wave functions in the center of 
the first Brillouin zone. We find that the strongest in¬ 
terference patterns are not obtained when the excitation 
wavelength is equal to the nearest neighbor separation a^. 
Instead, the symmetry of the lattice favors wavelengths 
that are 13.4 % shorter than ao- The discrepancy hints 
at the non-trivial symmetry properties of the honeycomb 
lattice and is explained by means of solid state theory, 
testifying to the applicability of our ripple tank for edu¬ 
cational purposes. 

The paper is organized as follows: section [I| describes 
the technical details of the setup, in section |n| and III 
we measure its physical properties, and in section IV 
the favorable excitation wavelengths are extracted from 
the discrete Fourier transform of a finite honeycomb lat¬ 
tice. Section [V| targets the correspondence between water 
waves in the ripple tank and electronic wave functions in 
graphene. We close with summary and outlook. 


I. SETUP 

Figure gives an overview of the modified ripple tank 
setup. A loudspeaker is used as a precise and inexpen¬ 
sive source of excitation.Sl It is mounted face-down inside 
a plexiglass hood. An aluminium rack holds both the 
hood and an interchangeable source plate 295 x 295 x 
4.2 mm^ in size. The source plate seen in Fig. [^is pat¬ 
terned with a honeycomb lattice of holes. Directly under¬ 
neath the source plate we place a commercially available 
ripple tank (Leybold, No. 401 501, with see-through wa¬ 
ter trough, tilted mirror and observation screen). The 
distance between source plate and water surface is only a 



FIG. 1. (a): Photograph of the observation screen picturing 

the ripple tank’s surface while driving air through a plexiglass 
plate with 660 hexagonally arranged holes. The next neighbor 
distance is ao = 1 cm. The excitation frequency that generates 
the sharpest patterns corresponds to a wavelength in water 
of A = V3/2cm. (b): like in (a), but with an additional 

phase difference between stroboscope and excitation of tt. The 
centers of the dark hexagons correspond to the bright spots 
in (a). Over time, the wave oscillates back and forth between 
these patterns, (c/d): Simulation of the 2D wave pattern 
obtained by solving the wave equation for 660 time-harmonic 
sources in a water-like environment, see Eq. Parameters 
like in (a) and (b), respectively. The white lines in (c) indicate 
the honeycomb pattern of the source plate. 


few mm. A white LED mounted in front of the center of 
the loudspeaker illuminates the ripple tank surface and 
can be used as stroboscope to observe a wave pattern at 
fixed phase. To that end, the LED is typically powered 
for 1 % of time of each wave cycle, thus allowing the 
observation of the wave during a 4° phase interval. 

Both the stroboscope’s and the loudspeaker’s frequen¬ 
cies, relative phase-shift, wave form and amplitude are 
controlled with a LabView program on a standard PC 
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FIG. 2. Setup of the ripple tank showing a 2D honeycomb 
wave pattern. Sound waves travel from the loudspeaker (top) 
through the hexagonally patterned source plate onto the rip¬ 
ple tank’s surface. The image is created by a stroboscopic 
LED (top) that projects the wave pattern via the mirror onto 
the screen. 


via an USB Audio Adapter (LogiLink, UA0078). Be¬ 
tween PC and loudspeaker or stroboscope, amplifiers are 
used to scale the signals appropriately. The loudspeaker 
requires between 0.3 and 3V AC input voltage, the LED 
4V DC. The loudspeaker receives a sinusoidal input sig¬ 
nal, the stroboscope’s waveform is set to rectangular. 

In combination with the loudspeaker, the hole pattern 
of the source plate thus creates a controlled spatial pat¬ 
tern of point sources of sinusoidal-in-time airblasts di¬ 
rected downward onto the surface of the water.l^ Since 
ripple crests/valleys act on the illuminating light as con¬ 
verging/diverging lenses,!^ the resulting wave pattern on 
the water surface becomes visible as bright and dark ar¬ 
eas projected onto the observation screen. 


II. SINGLE SOURCE 


We performed a series of single source experiments to 
determine the dispersion relation of water waves in our 
setup. The water height in the tank was kept constant at 
h = l.S cm. The wavelength lies between 0.5 cm and 2.5 
cm. Since wavelength and water height are of the same 
order of magnitude, the dispersion relation for ripples is 
to a good approximatiorP given by 

= gk+ - (1) 

p 



EIG. 3. (a): Photograph of the observation screen for a 

single point source plate at A = 1 cm. (b): Measured (circles) 
and calculated (line) dispersion relation in the ripple tank, 
(c): Photograph of the observation screen for a double point 
source plate at A = 1 cm. (d): Simulation of the 2D wave 
pattern created by two time-harmonic sources in a water-like 
environment. 


where uj is the angular frequency of the ripple, k = 27r/A 
the wavevector, p the density of the water (Ig/cm^), S 
the surface tension (73g/s^), and g the absolute value of 
the gravitational acceleration (981cm/s^). 

For single source experiments, we fabricated a source 
plate with a single hole 2 mm in diameter. The frequency 
was swept from 5 to 40 Hz. With a digital camera, we 
took pictures of the observation screen. The radial dif¬ 
ference between two adjacent bright circles in Fig. [^a) 
corresponds to the wavelength. Fig.j^b) shows our mea¬ 
sured wavelengths in comparison to Eq. For frequen¬ 
cies above 10 Hz, we find a good agreement. The devia¬ 
tion at low excitation frequencies could be caused by op¬ 
erating the loudspeaker outside the specified range which 
might introduce nonlinearities in the excitation. Fitting 
the data using least mean squares, we extract a mea¬ 
surement uncertainty of cr =1 Hz for each data point, 
assuming that all uncertainties are statistical in nature. 
Measurement uncertainties comprise calibration errors of 
the gauge scale on the observation screen, reading errors 
of the wavelength and calibration errors of the frequency 
generator. Thus, for example at 30 Hz, the wavelength 
of the water waves generated by our setup is known with 
0.02 cm precision. 
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III. TWO SOURCES 

Using a plate with two holes, we can emulate a double 
slit experiment. We choose this well-known geometry to 
compare measurement with simulation as well as to study 
the effect the finite hole diameter has on the pattern. 

The aisles of constructive and destructive 
interference can be read from Fig. |^c). Starting 
with the solution to the wave equation of a 2D 
time-harmonic source for a radially propagating wave 

^ir,k,co,t) = -^{H^[kry-^ + Hl[-kr]e-^^^) (2) 

with Hq the zero order Hankel function of the second 
kind, k the absolute value of the wave vector, r the radial 
distance from the source, uj the radial frequency, and t 
the time, we performed a Mathematica-based simulation 
of the 2D wave pattern. The result of these calculations 
is shown in Fig. |^d). 

To estimate the effect of the finite hole size, we compare 
a simulated wave function 'ipi{x,y) with two point-like 
sources to a wave function where each source 

is replaced by five sources at the four endpoints and the 
center of a crosshair whose diameter corresponds to the 
diameter of the holes. As a figure of merit, we introduce 
the correlation between ' 0 i(x,^) and 2 p 2 {x,y) given by: 

^ ^ ^2{x,y)) 

Using a distance between the sources of ao = 1 cm and a 
hole diameter of 2 mm we find C to be 99.7 %. Since these 
parameters are also convenient to realize experimentally, 
we use them with confidence in all subsequent source 
plates. 


IV. HONEYCOMB LATTICE 

The symmetry properties of the honeycomb lattice can 
be inferred by performing the Fourier transform of the 
wave pattern. We exploit that this transform can be 
expressed as the product of the Fourier transform of a 
single time harmonic point source and that of the lattice. 
We perform a discrete Fourier transform of a finite lattice 
and compare the results to key parameters known from 
the Fourier transform of an infinite lattice, which are 
deduced from the pseudopotential method for a crystal 
with a two-atom basis 

We created a 42 point lattice with uq = 1 cm for the 
discrete transform. Figure [^a) pictures this lattice in 
real space. Figure |^b) shows the squared modulus of its 
discrete Fourier transform. As expected from the pseu¬ 
dopotential method we find a pattern with two classes of 
points with an intensity ratio of 1 to 4. Purple arrows 



FIG. 4. (a): Honeycomb pattern with 42 holes: nearest- 

neighbor distance ao is 1 cm. Dashed black lines indicate the 
unit cell of the lattice, (b): 2D discrete Fourier transform 
of the pattern depicted in a. Dark spots correspond to high 
amplitude. Purple arrows indicate the six lowest reciprocal 
lattice vectors, the black arrow indicates the next highest re¬ 
ciprocal lattice vector, (c): 2D discrete Fourier transform of 
a single time harmonic source with wavelength A = V3/2cm. 
(d): 2D discrete Fourier transform of 42 time harmonic point 
sources with wavelengths of A = v^/2cm. 

indicate the 6 shortest reciprocal lattice vectors, 6 i, 62 , 
— 61 , — 62 , — 1)2 and —bi -h ^ 2 , with 

and correspond to the weak intensity points in the pat¬ 
tern. The black arrow indicates the next longer recipro¬ 
cal lattice vector 61+62 . The radii of the dashed cir¬ 
cles in pink and black correspond to the /c-vector lengths 
favoured by the pattern, while all other lengths are sup¬ 
pressed. The finite lattice size employed in the calcula¬ 
tion manifests itself in the width of the points. 

The discrete 2D Fourier transform of a single time har¬ 
monic point source is shown in Fig. [^c) for the spe¬ 
cific wavelength of v^/ 2 cm. We see a circle of radius 
47r/v^cm“^. 

The final result, the 2D Fourier transform of the 42 
time harmonic point sources is shown in Fig. [^d), which 
pictures the product of the Fourier transforms shown in 
Fig. [^b) and Fig. [^c). 

Consequently, the discrete 2 D Fourier transform of 
Fig. B a) which was taken at an excitation wavelength 
of A = \/3/2cm yields six points on a circle in k-space. 
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analogous to Fig. [^d). 

This analysis illustrates that the resulting wave pat¬ 
tern in real space is, at the given wave length, essentially 
a superposition of six plane wave states. In general, the 
excitation wavelength, which determines the radius of 
the circle in Fig. |^c), therefore serves as a tool to se¬ 
lect the finite set of plane wave states that superimpose 
in real space. Most importantly, it is not the correspon¬ 
dence between the nearest neighbor separation and the 
wave length in real space that gives rise to strong inter¬ 
ference patterns, as one could naively think. Instead, the 
wavelength of strongest interference A = \/3/2ao corre¬ 
sponds to half the side length of the rhombus as indi¬ 
cated by the black lines in Fig. [^a). Thus, the lattice 
ripple tank experiment leads to the conclusion that the 
rhombus and not the honeycomb is the unit cell of the 
graphene latticc.^^ The rhombus is the smaller and thus 
more fundamental unit that any honeycomb lattice can 
be constructed from. 


V. CORRESPONDENCE TO GRAPHENE 

In the present experimental setting, all sources are in 
phase, since they are excited by a common air pulse. This 
means that as far as the correspondence to graphene is 
concerned, we are limited to mimicking wave functions 
consisting of counterpropagating waves at the F- point. 



1st Brillouin zone 


FIG. 5. Cut through the 2D bandstructure of graphene 
along the kx direction (schematic drawing). On the x-axis, 
the symmetry points and the first Brillouin zone are indicated. 
The dashed lines indicate the lifting of the degeneracy due to 
the atomic potential. 


The strongest resonance in the ripple tank is_ excited 

it cor- 
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with A = \^/2ao. As discussed in section 
responds to k -vectors of length k = . There is a 

superposition of six basis vectors of the reciprocal lat¬ 
tice with this k-vector length, see Fig. Bb). Labeling 
them consecutively by Ai through Aq in clockwise direc¬ 
tion and solving the Schrodinger equation for graphene 
in the emty lattice approximation, we can draw the 2D 
bandstructure that consists of six parabolic dispersison 
relations all crossing at the F- pointFigure shows a 
cut along the A 2 - A 5 direction which coincides with the 
kx direction. Beyond the empty lattice approximation. 


the six-fold degeneracy at the F- point is lifted by tak¬ 
ing the atomic potential into account. The lowest-lying 
band, which is the symmetric combination of all six basis 
states with wave-vectors Ai to Ae and therefore the only 
combination without phase differences between the waves 
corresponds to the symmetry of the water waves and will 
therefore exhibit the same symmetry as the honeycomb 
pattern shown in Fig. 


VI. SUMMARY AND OUTLOOK 

With a modification that provides periodic sources, we 
have extended the classic ripple tank experiment to il¬ 
lustrate symmetry properties of the honeycomb lattice. 
The experimentally observed patterns can be reproduced 
closely by numerical simulations. A next step could be 
to introduce phase shifts between the individual point 
sources. However, this requires a completely new exci¬ 
tation technique. A more promising application for the 
present technique is to mimic electron waves in cavities 
In this case, the apparatus could not only be used to make 
current science accessible to students or the general pub¬ 
lic, but might also prove to be a useful modelling tool 
for microstructured geometries that are costly and time- 
consuming to fabricate. 
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